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Introduction 

'he  proposed  fuel,  RJ-5,  is  a blend  of  several  fused  ring  hydro- 
carbons of  low  carbon/hydrogen  ratio.  This  property  imports  high 
density  and  high  volumetric  heat  of  combustion  to  the  fuel,  features  which 
make  it  a promising  candidate  for  volume  - limited  applications  such  as 
the  ASALM  program.  Unfortunately,  the  need  for  improvements  in 
several  of  its  physical  properties,  such  as  freezing  point  and  viscosity, 
as  well  as  production -related  difficulties,  have  necessitated  modification 
of  the  fuel's  isomer  ratio  and  the  use  of  additives. 

It  is  well  established  that  small  quantities  of  additives  can  have 
major  effects  on  combustion  parameters  such  as  flame  velocity'  , first 
and  second  ignition  delays'^  , and  hot-wire  ignition  temperatures^. 

Many  of  these  effects  are  difficult  to  predict  a prio ri  and  to  explain 
a posteriori.  For  this  reason,  considerable  interest  has  been  generated 
in  the  development  of  new  and  innovative  techniques  to  monitor  the  effect 
of  variation  in  fuel  composition  and  to  obtain  a better  understanding  of 
fuel  combustion  kinetics  by  using  additives  to  tailor  fuel  performance  to 
engine  needs. 

One  promising  tool  for  combustion  research  is  the  chemical  shock 
tube.  With  this  method,  fuel-air  mixtures  can  be  rapidly  heated  to  a 


preselected  temperature  and  pressure,  and  their  combustion  can  be 
studied  free  from  the  complicating  temperature,  pressure,  and  density 
gradients  associated  with  flame  studies. 


Two  previous  shock  tube  studies  on  the  combustion  kinetics  of 
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RJ-5  have  been  carried  out  by  Shell  , and  Martin  Marietta  Laboratories 
(MML)^'  to  measure  ignition  delays,  t . The  goals  of  these  preliminary 
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studies  were  to:  (1)  investigate  the  effect  of  reaction  conditions  on  r, 
and  (2)  determine  the  influence  of  viscosity- reducing  additives.  While 
extensive  data  on  the  variation  of  r were  obtained  with  the  experimental 
techniques  employed,  it  is  not  clear  how  these  data  would  correlate 
with  engine  test  firing  results  because  r is  merely  an  indicator  of  the 
onset  of  combustion  and  is  not  directly  related  to  combustion  duration 
or  efficiency. 

Because  of  the  uncertainties  associated  with  ignition-delay  measure- 
ments, a more  extensive  study  was  carried  out  at  MML  to  establish 
procedures  for  monitoring  the  entire  combustion  profile  for  RJ-5  vapor  -- 
air  mixtures  rather  than  just  ignition^.  In  this  study,  the  infrared 
emissions  emanating  from  a hot  shock -initiated  combustion  zone  were 
filtered  to  yield  CO^  and  H^O  time  profiles.  From  these  data,  the  effect 
of  variation  in  the  reaction  temperature,  pressure,  and  fuel/air  ratio 
on  the  combustion  rate  and  duration  were  obtained. 

It  was  recognized,  however,  that  combustion  in  real  systems  is 

heterogeneous,  i.  e.  , ignition  of  fuel  vapor  occurs  first,  followed  by  the 

much  slower  burning  of  a fine  mist  of  liquid  droplets  in  the  hot  combustion 

zone.  This  aspect  of  the  combustion  problem  is  more  significant  for 

RJ-5  than  for  JP-4  or  7 because  RJ-5's  higher  viscosity  leads  to  much 

(7) 

larger  droplets  due  to  air  blast  shearing  in  rocket  engines'  . 

The  present  work  has  been  undertaken  to  study  the  heterogeneous 
combustion  of  RJ-5  air  mist.  To  study  this  complex  but  technically 
more  realistic  situation,  we  have  added  an  ultrasonic  nebulizer  to  the 
downstream  end  of  a shock  tube  (Fig.  1).  With  this  device,  a fine  mist 
of  controlled  particle  size  and  concentration  is  generated  near  the  shock 
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32  I.  Driver  Gas  Tank  Spring  Loaded  Knife 


General  schematic  of  nebulizer -coupled  shock  tube 


tube  end  wall.  A shock  wave  is  initiated  through  air,  compressing  and 
heating  the  mist  and  causing  combustion  of  the  atomized  fuel.  The  con- 
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struction  details  of  the  nebulizer-coupled  shock  tube  system  have  already 
been  discussed 

The  present  report  will  be  concerned  with  the  formation,  stability, 
and  characterization  of  the  precombustion  shock-induced  RJ-5  air  mis1. 
Shock  Wave  Production 

Since  combustion  in  heterogeneous  systems  must  be  preceded  by 
vaporization,  it  is  inherently  slower  than  simple  vapor  phase  combustion. 
This  difference  would  probably  lead  to  combustion  times  that  would 
exceed  the  reflected  shock  wave  duration  of  the  existing  MML  system. 

To  a first  approximation,  the  useable  duration  of  the  reflected 
shock  zone  is  proportional  to  the  shock  tube  length.  However,  wall  drag 
(in  small  diameter  tubes)  and  inadequate  diaphragm  openings  are  impor- 
tant parameters  in  causing  reflected  shock  decay.  It  was  felt  that  these 
deficiencies  could  not  be  overcome  with  the  existing  shock  tube.  Conse- 
quently, a new  tube  has  been  constructed  to  eliminate  thest  two  difficulties. 

To  ensure  uniform  and  rapid  diaphragm  rupture,  a spring-loaded, 
double -bladed  knife  has  been  installed  into  the  driver  sections.  Dia- 
phragm rupture  is  initiated  by  release  of  the  cocked  knife  which  accel- 
erates downstream  and  cuts  the  diaphragm  in  a quadrant  petaling  pattern. 
The  knife  is  rapidly  retracted  from  the  foil  region  by  spring  tension, 
producing  an  unhindered  path  for  driver  gas  expansion,  and  hence,  a 
longer  lifetime  for  the  reflected  zone.  Aluminum  diaphragms  ratain 
their  post-shock  configuration  and,  readily  exhibit  the  quadrant  pattern 
(Fig.  2). 
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To  minimize  the  effect  of  wall  drag,  1.  75  in.  i.  d.  square  polished 
aluminum  pipe  was  used  for  the  shock  tube  construction  instead  of  the 
1 in.  i.  d.  circular  tubing  previously  employed.  This  corresponds  to  a 
fourfold  increase  in  shock  tube  cross-sectional  area,  which  should 
s ignificantly  reduce  drag  effects.  Furthermore,  care  was  taken  to 
ensure  that  shock  tube  attachments,  e.  g.  , windows,  valves,  pressure 
transducers,  etc.  , were  flush-mounted  to  limit  turbulent  effects. 

Figure  3 shows  clearly  that  the  reflected  shock  zone  is  stable 
beyond  1500  p sec.  As  will  be  demonstrated  later,  this  time  interval 
should  be  ample  for  vaporization  and  combustion  of  the  fuel  mist. 
Aerosol  Characte rization 

The  coupling  of  an  ultrasonic  nebuli  and  shock  tube  has  not 
previously  been  attempted.  For  this  reason,  cedures  for  mist 

characterization  in  shock  tubes  have  not  been  fu1  } established.  Hence, 
the  factors  affecting  aerosol  properties  are  discussed  in  some  detail 
below. 

1 ) Particle  Concentration 

Upon  initiation  of  nebulization  in  the  preshocked  combustion  zone, 
a fine  air  mist  will  be  formed.  However,  droplet  settling  will  begin 
immediately.  This  process  is  depicted  in  the  figure  shown  below. 


SHOCK  TUBE  PRESSURE 


at  2000°K 


(b)  Line  drawing  of  above  experiment 

a)  Incident  shock  arrival  at  observation  station 

b)  Reflected  shock  arrival  at  same  station. 


apidly  be  achieved;  thus,  the  mist  concentration  will  be  independent 


e expression 


(droplets) 


where  and  kg  are  the  nebulization  and  settling  rate  constants, 
respectively.  Thus,  the  steady  droplet  concentration,  (droplets),  is 


(droplets) 


The  value  of  k~.  can  be  obtained  from  the  time  required  to  nebulize  a 


known  quantity  of  fuel  under  a flowing  condition  where  the  droplets 


cannot  return  to  the  transducer  surface.  The  settling  rate  constant 


can  be  calculated  as  discussed  below.  Thus,  both  factors  affecting  the 


steady  state  particle  concentration  can  be  readily  obtained 


The  settling  velocity  of  spherical  particles  in  a stagnant  air  mass 

(9) 

is  given  by  the  "tranquil  settling  expression"'  . 


where  d is  the  droplet  diameter,  77  is  the  viscosity  of  the  medium,  and 


p and  p are  the  densities  of  the  particles  and  settling  medium,  re 


spectively.  For  RJ-5  droplets  in  air,  p . «p 


be  neglected.  Also,  the  viscosity  of  air  is  essentially  independent  of 


pressure  over  the  region  of  interest.  The  calculated  settling  rates  of 


the  droplet  sizes  of  interest  to  the  present  study,  those  of  diameter 


0.5  - 50  p,  are  given  in  Table  I 


Of  significance  here  is  the  fact  that  the  particle  settling  time, 
indicated  1/C,  is  a strong  function  of  particle  diameter.  This  has  two 


Settling  Rates  of  RJ-5  Droplets  in  Air 


• It  will  be  necessary  to  vary  the  irradiation  power  to  maintain 


the  same  droplet  concentration  as  the  mean  droplet  diameter 


i3  varied  via  the  irradiation  frequency;  and 
• The  steady-state  particle  distribution  of  the  mist  will  differ 
somewhat  from  that  of  the  droplets  as  they  leave  the  transducer 


surface.  Fortunately,  ultrasonic  nebulizers  have  fairly  narrow 


distribution  curves.  Hence,  this  difference  should  not  pose  any 


oblems  to  mist  characterization 


Particle  Size  Distribution 


In  addition  to  the  ease  with  which  it  could  be  coupled  to  the  shock 


tube,  an  important  reason  for  the  choice  of  the  ultrasonic  nebulizer  for 


producing  the  fuel  aerosol  was  its  ability  to  vary  the  mean  particle  size 


1 


Furthermore,  as  mentioned  previously,  the  distribution  is  relatively 
narrow.  This  is  exemplified  in  the  data  given  by  Dentor/^'  ibelow). 


3MHz 


1 MHz 


Droplet  Size  in  Microns 


However,  the  theory  of  ultrasonic  nebulizers  has  not  been  adequately 
developed  to  account  for  changes  in  liquid  physical  properties,  such  as 
viscosity,  surface  tension,  etc.  Consequently,  it  will  be  necessary  to 
experimentally  determine  the  droplet  size  distribution  curve  for  RJ-5 
as  a function  of  the  irradiation  frequency.  Of  the  various  methods 
available,  the  approach  taken  by  Kousaka'  seems  ideally  suited  for 
nebulizer  characterization.  In  this  method,  particle  size  distribution 
is  determined  by  the  time  and  number  of  particles  settling  through  the 
depth  of  focus  of  an  ultramicroscope.  This  method  has  been  successfully 
employed  to  characterize  water  aerosols  in  the  range  oi  0.  5 to  10  p. 

The  equipment  necessary  to  carry  out  this  procedure  is  presently 
available  at  MML. 

3 ) Particle  Vaporization 

Prior  to  shock  arrival  in  the  mist  region,  the  RJ-5  liquid -vapor 
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equilibrium  will  already  be  established.  Consequently,  mass  loss 
through  vaporization  will  not  be  a significant  problem  in  the  preshock 
mist.  However,  vaporization  will  begin  immediately  upon  shock  heating. 
The  total  time  required  to  vaporize  the  droplets  in  a stagnant,  non- 
combustion atmosphere,  e.  g.  , a shock -heated  nitrogen,  rather  than 
air-mist,  can  be  readily  calculated  from  the  relationship  (refer  to 
Appendix  I) 

r = P H r 2/(2 A f T - T,1 ) 
v v o ' L m dj 

where  Hv  is  the  heat  of  vaporization,  rQ  the  initial  drop  radius,  A is  the 

thermal  conductivity  of  the  atmosphere,  and  T and  T , are  the  air  and 

m a 

drop  temperatures,  respectively.  For  instance,  a 50  p (dia. ) drop  of 
RJ-5,  initially  at  300°K  and  shock  heated  to  1273°K  (1000°C),  will 
completely  vaporize  in  about  1400  p sec.  , which  is  about  the  lifetime 
of  the  reflected  shock  zone  (refer  to  Fig.  3).  However,  when  surrounded 
by  flame  in  a combustion  environment,  the  vaporization  time  will  be 
substantially  less.  Therefore,  it  is  unlikely  that  combustion  of  50  p 
RJ-5  droplets  would  be  incomplete  in  the  observation  time  available. 

4)  Axial  Stability 

Up  on  arrival  of  the  incident  shock  wave  in  the  mist  region,  the 
aerosol  particles  will  suddenly  come  under  the  influence  of  a frictional 
force  or  drag,  which,  given  time,  will  accelerate  them  to  supersonic 
shock  speed.  Fortunately,  however,  the  accelerating  force  is  instantly 
neutralized  when  the  particles  collide  with  the  stagnant  reflected  shock 
zone.  Although  the  time  between  the  incident  and  reflected  shock 
arrival,  t , is  short  --  about  50  p sec.  (refer  to  Fig.  3),  the  particles 
will  have  a finite  velocity  in  the  direction  of  the  end  plate  when  they 
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enter  the  reflected  :.onr.  Similarly,  they  will  then  decelerate  dvie  to 
frictional  drag  with  the  stationary  reflected  shock  gasses.  Since  it  is 


important  that  the  burning  RJ-5  mist  remain  in  the  viewing  region 
(located  about  1 cm  from  the  end  plate)  and  not  be  swept  into  it  during 
the  observation  time,  an  analysis  of  the  mist  axial  drift  is  presented 
below. 

According  to  Stokes,  a spherical  particle  moving  in  a stationary 
gas  stream,  or  conversely,  a gas  stream  flowing  over  a stationary 
particle,  will  exert  a frictional  force  on  the  particle  described  by  the 
relationships  ^ ^ 


Force  = 6»nrv 


-m 


dv 

dt 


where  r and  m are  the  particle  radius  and  mass,  respectively,  and  v 
is  the  particle  velocity  relative  to  the  gas.  This  relationship  leads  to 
the  expression  for  the  particle  speed,  s,  relative  to  the  side  viewing 
port  (Appendix  II).  . - ~ 


where  p is  the  droplet  density,  and  "a"  is  the  incident  shock  velocity. 

The  velocity,  s,  of  a 50  p dia.  KJ  -5  droplet  exposed  to  a Mach  2 incident 
shock  wave  would  be  1,180  cm/sec.  after  50  p sec.  Thus,  during  this 
50  p sec.  , it  would  have  moved  only  0.  35  mm  from  its  original  preshock 
position.  Furthermore,  the  reflected  shock  gas  viscosity,  and  hence 
its  decelerating  force,  would  be  substantially  higher  than  in  the  incident 
zone  dvie  to  the  higher  temperature.  Consequently,  droplets  of  approxi- 
mately this  size  should  remain  in  the  observation  region  during  combustion. 
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The  rate  of  heat  transfer  to  a spherical  drop  of  RJ-5  from  the 

( 1 2) 

surrounding  atmosphere,  dg/dt,  is' 
dg/dt  = 4 r A(Tm  - Td) 

where  X is  the  thermal  conductivity  of  the  atmosphere.  It  is  is  assumed 
that  all  the  heat  transferred  is  used  to  vaporize  RJ-5  at  constant  temp- 
erature T^,  then 

dg/dt  = -Hv  dm/dt  = 4 7rpr^  Hy  dr/dt,  where 
where  is  the  heat  of  vaporization  per  unit  weight,  and  p is  the 
density  of  RJ-5.  Combining  these  two  equations  leads  to 
dr/dt  = - X (Ta  - Td)/(  P r Hy) 

which,  upon  integration,  leads  to  an  expression  for  the  total  vaporization 
time, 

'v  ■ "Hvro2/[|2i(Ta  -Td>]- 


